TRUNCATED DECOMPOSITIONS AND FILTERING METHODS 
WITH REFLECTIVE/ ANTI-REFLECTIVE BOUNDARY 
CONDITIONS: A COMPARISON 

C. TABLING POSSIO* 

Abstract. The paper analyzes and compares some spectral filtering methods as truncated 
singular/eigen- value decompositions and Tikhonov/Re-blurring regularizations in the case of the re- 
cently proposed Reflective 1181 and Anti-Reflective |21| boundary conditions. We give numerical 
evidence to the fact that spectral decompositions (SDs) provide a good image restoration quality 
and this is true in particular for the Anti-Reflective SD, despite the loss of orthogonality in the 
associated transform. The related computational cost is comparable with previously known spectral 
decompositions, and results substantially lower than the singular value decomposition. The model 
extension to the cross-channel blurring phenomenon of color images is also considered and the related 
spectral filtering methods are suitably adapted. 
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1. INTRODUCTION. In this paper we deal with the classical image restora- 
tion problem of blurred and noisy images in the case of a space invariant blurring. 
Under such assumption the image formation process is modelled according to the 
following integral equation with space invariant kernel 

g{x)= h{x-x)f{x)dS: + Tj{x), xeR^, (1.1) 



where / denotes the true physical object to be restored, g is the recorded blurred 
and noisy image, rj takes into account unknown errors in the collected data, e.g. 
measurement errors and noise. 



As customary, we consider the discretization of (1.1 1 by means of a standard 2D 
generalization of the rectangle quadrature formula on an equispaced grid, ordered 
row-wise from the top-left corner to the bottom-right one. Hence, we obtain the 
relations 

5, = ^ /i,_j/,+?7„ ^eZ^ (1.2) 

in which an infinite and a shift-invariant matrix A^c — [^j-j](i,j)=((ii.i2).(ji.i2))' ^ 
two-level Toeplitz matrix, is involved. 



In principle, (1.2 1 presents an infinite summation since the true image scene does 
not have a finite boundary. Nevertheless, the data gi are clearly collected only at a 
finite number of values, so representing only a finite region of such an infinite scene. In 
addition, the blurring operator typically shows a finite support, so that it is completely 
described by a Point Spread Function (PSF) mask such as 



l-PSF 



\hi,^^■ (1.3) 
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where hi-^^i^ > for any ^1,12 and J2i=-q^i = 1; * — (*ii*2), 9 — {<li,q2) (normaliza- 
tion according to a suitable conservation law). 



Therefore, relations (1.2 1 imply 



q 

^hsfi-s + Vt^ zi = l,...,ni,i2 = l,...,n2, (1.4) 

s=-q 



where the range of collected data defines the so called Field of View (FOV) . 



Once again, we are assuming that all the involved data in (1.5), similarly to (1.2 1, are 



reshaped in a row-wise ordering. In such a way we obtain the linear system 

Af^g-V (1-5) 

where A S K^(")x^("+29) jg ^ finite principal sub-matrix of Aoo, with main diagonal 
containing ho Q, f g I|J^("+29)^ 9,V & and with N{m) = mim2, for any two- 

index m = (toi, ?Ti2). 

Such a reshape is considered just to perform the theoretical analysis, since all the 
deblurring/denoising methods are able to deal directly with data in matrix form. 



For instance, it is evident that the blurring process in (1.4) consists in a discrete 
convolution between the PSF mask, after a rotation of 180°, and the proper true 
image data in 



■(7i + l,...,Tii+(ji,i2 = — g2 + l,---,n2+g2 



Hereafter, with a two-index notation, we denote hy F — \ fi, iJ . , , the 

' ' ^ ^1 1 ^2 J 4^ — 1^. . . .m ,22 — 1,. . . ,n2 

true image inside the FOV and by G = [5ii,i2]ij=i m 1-2=1 recorded image. 



Thus, assuming the knowledge of PSF mask in ( 1.3 1 and of some statistical properties 



of 77, the deblurring problem is defined as to restore, as best as possible, the true 



image F on the basis of the recorded image G. As evident from (1.4), the problem is 
undetermined since the number of unknowns involved in the convolution exceeds the 
number of recorded data. Boundary conditions (BCs) are introduced to artificially 
describe the scene outside the FOV: the values of unknowns outside the FOV are fixed 
or are defined as linear combinations of the unknowns inside the FOV, the target being 



to reduce (1.5) into a square linear system 

AJ^g-T] (1.6) 
with A^ e ^Nin)xN(n)^ (ni,n2), N{n) = nin2 and £ 

The choice of the BCs does not affect the global spectral behavior of the matrix. 
However, it may have a valuable impact both with respect to the accuracy of the 
restored image and to the computational costs for recovering / from the blurred 
datum, with or without noise. 

Notice also that, typically, the matrix A is very ill-conditioned and there is a significant 
intersection between the subspace related to small eigen/singular values and the high 
frequency subspace. Such a feature requires the use of suitable regularization methods 
that allow to properly restore the image F with controlled noise levels [121 1131 El El] , 
among which we can cite truncated SVD, Tikhonov, and total variation [12 1 IT4 l [24], 
Hereafter, we focus our attention on special case of PSFs satisfying a strong symmetry 
property, i.e., such that 



= hi for any i = —5, ...,(?. 



(1.7) 
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This assumption is fulfilled in the majority of models in real optical applications. 
For instance, in most 2D astronomical imaging with optical lens [5 the model of 
the PSF is circularly symmetric, and hence, strongly symmetric; in the multi-image 
deconvolution of some recent interferometric telescopes, the PSF is strongly symmetric 
too [6] . Moreover, in real applications when the PSF is obtained by measurements (like 
a guide star in astronomy) , the influence of noise leads to a numerically nonsymmetric 
PSF, also when the kernel of the PSF is strongly (or centro) symmetric. In such a case, 
by employing a symmetrized version of the measured PSF, comparable restorations 
are observed [TSIIT]. 

The paper is organized as follows. In Section [2] we focus on two recently pro- 
posed BCs, i.e., the Reflective [18^ and Anti-Reflective BCs [21] and their relevant 
properties. Section |3] summarizes some classical filtering techniques as the truncated 
singular/eigen- values decomposition and the Tikhonov method. The Re-blurring 
method [HllS] is considered in the case of Anti-Reflective BCs and its re-interpretation 
in the framework of the classical Tikhonov regularization is given. In Section [4] the 
model is generalized for taking into account the cross-channel blurring phenomenon 
and the previous filtering methods are suitable adapted. Lastly, Section [5] deals with 
some computational issues and reports several numerical tests, the aim being to com- 
pare the quoted filtering methods and the two type of BCs, both in the case of gray- 
scale and color images. In Section [6] some conclusions and remarks end the paper. 

2. BOUNDARY CONDITIONS. In this section we summarize the relevant 
properties of two recently proposed type of BCs, i.e., the Reflective [18] and Anti- 
Reflective BCs [H] . Special attention is given to the structural and spectral properties 
of the arising matrices. In fact, though the choice of the BCs does not affect the global 
spectral behavior of the matrix A, it can have a valuable impact with respect both 
to the accuracy of the restoration (especially close to the boundaries where ringing 
effects can appear), and the computational costs for recovering the image from the 
blurred one, with or without noise. 

Moreover, tanking into account the scale of the problem, the regularization methods 
analysis can be greatly simplifled whenever a spectral (or singular value) decomposi- 
tion of A is easily available. This means that the target is to obtain the best possible 
approximation properties, keeping unaltered the fact that the arising matrix shows 
an exploitable structure. For instance, the use of periodic BCs enforces a circulant 
structure, so that the spectral decomposition can be computed efficiently with the 
fast Fourier transform (FFT) [8J. Despite these computational facilities, they give 
rise to signiflcant ringing effects when a significant discontinuity is introduced into 
the image. 

Hereafter, we focus on two recently proposed boundary conditions, that more care- 
fully describe the scene outside the FOV. 

Clearly, several other methods deal with this topic in the image processing literature, 
e.g. local mean value [22] or extrapolation techniques (see [T7] and references therein). 
Nevertheless, the penalty of their good approximation properties could lie in a linear 
algebra problem more difficult to cope with. 

2.1. REFLECTIVE BOUNDARY CONDITIONS. In [18, Ng et al. ana- 
lyze the use of Refiective BCs, both from model and linear algebra point of view. The 
improvement with respect to Periodic BCs is due to the preservation of the continuity 
of the image. In fact, the scene outside the FOV is assumed to be a refiection of 
the scene inside the FOV. For example, with a boundary at xi = and 2:2 = the 
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reflective condition is given by f{±xi,±X2) = f{xi,X2)- 
More precisely, along the borders, the BCs impose 

fii,l-i2 = /ii.i2i fii,n2+i2 — /ji,n2 + l-j2' fol' *1 = 1; ■ • ■ ' '^l' *2 = 1, ■ • ■ , 92 

fi-ii,i2 — fni+ii,i2 = /ni+i-ii,j2i for any ii = 1, . . . ^2 = 1, • ■ • ,"2, 

and, at the corners, the BCs impose for any ii = 1, . . . , gi, Z2 = 1, ■ ■ ■ , 92 



fl — 11,1 — 12 /ii,i2' 
fl — ii,n2+i2 ^ /ii,n2 + l — i2i 



fni+ii,n2+i2 /rii +1 — ii ,n2 + l— 12 ' 

fni+ii,l — i2 /rii+1 — 11,12 > 



i.e., a double reflection, first with respect to one axis and after with respect to the 
other, no matter about the order. ^ 

As a consequence the rectangular matrix A is reduced to a square Toeplitz-plus- 
Hankel block matrix with Toeplitz-plus-Hankel blocks, i.e., A„ shows the two-level 
Toeplitz-plus-Hankel structure. Moreover, if the blurring operator satisfies the strong 
symmetry condition ( |1.7[ ) then the matrix An belongs to DCT-III matrix algebra. 
Therefore, its spectral decomposition can be computed very efficiently using the fast 
discrete cosine transform (DCT-III) [SS]. 

More in detail, let C„ = {A„ e M^(")><^("), n = {ni,n2),N{n) = ni?i2 | A„ = 
i?„A„i?^} be the two-level DCT-III matrix algebra, i.e., the algebra of matrices that 
are simultaneously diagonalized by the orthogonal transform 



R„ 



R„ 



5t,l ^^^J{s-l){t~l/2)TT 



(2.1) 



with 5s.t denoting the Kronecker symbol. 

Thus, the explicit structure of the matrix is An = Toeplitz(l/)+Hankel(cr(F), Ja{V)), 
with y = [Vo ^1 ... . . . 0] and where each 1^^, ii — 1, . . . , gi is the unilevel 
DCT-III matrix associated to the i^^ row of the PSF mask, i.e., Vi^ = Toeplitz(wij) -f 



Hankel(cr(wiJ, Ja(wiJ), with Vi^ = [^ii,o,-- 



, 0, . . . , 0]. Here, we denote by cr 



the shift operator such that a{vi-^) — . . . ,hi^ q^,Q, • • • ,0] and by J the usual flip 

matrix; at the block level the same operations are intended in block- wise sense. 
Beside this structural characterization, the spectral description is completely known. 



In fact, let / be the bivariate generating function associated to the PSF mask (1.3 1 
that is 



91 



92 



f{xi,X2) = + 2 ^ hs^fiCos{siXi) + 2 ^ /io,s2 cos(s2a;2) 

Sl=l S2 = l 

91 92 

+4 ^ ^ ^si,s2 cos(siXi)cos(s2a;2), 



(2.2) 



81 = 1 S2 = l 

then the eigenvalues of the corresponding matrix An G C„ are given by 

A«(^„) = /(x[7l,4'^^l), s=(si,S2), 4™'=^'^^;^, 

where Si = 1, . . . , ni, S2 = 1, . . . , n2, and where the two-index notation highlights the 
tensorial structure of the corresponding eigenvectors. 

Lastly, notice that standard operations like matrix- vector products, resolution of lin- 
ear systems and eigenvalues evaluations can be performed by means of FCT-III |18j 
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within 0(nin2 log(nin2)) arithmetic operations (ops). For example, by multiplying 
by d = [1, 0, . . . , 0]"^ both the sides of i?^A„ = A„i?^, it holds that 

[^"](il,i2) = [Rni^nel)](i,.^^■)/[Rnel](^,^^^), ii = 1, . . . , ni , is = 1 , . . . , ria , 

i.e., it is enough to consider an inverse FCT-III applied to the first column of An, 
with a computational cost of 0{nin2 log(nin2)) ops. 

2.2. ANTI-REFLECTIVE BOUNDARY CONDITIONS. More recently. 
Anti-reflective boundary conditions (AR-BCs) have been proposed in [21] and studied 
[21 [31 m ini [ini us • The improvement is due to the fact that not only the continuity of 
the image, but also of the normal derivative, are guaranteed at the boundary. This 
regularity, which is not shared with Dirichlet or periodic BCs, and only partially 
shared with reflective BCs, significantly reduces typical ringing artifacts. 
The key idea is simply to assume that the scene outside the FOV is the anti-reflection 
of the scene inside the FOV. For example, with a boundary at Xi = the anti-reflexive 
condition impose /(-a^i, a;2) - /(x^, 2:2) = -(/(xi, X2) - /(x^, 3:2)), for any X2, where 
x\ is the center of the one-dimensional anti-reflection, i.e., 

f{~xi,X2) = 2f{xl,X2) - /(a;i,a;2),for any X2- 

In order to preserve a tensorial structure, at the corners, a double anti-reflection, flrst 
with respect to one axis and after with respect to the other, is considered, so that the 
BCs impose 

f{-Xi,-X2) = 'ij{xl,xl) - 2f{xl,X2) - 2/(xi,X2) -I- /(a;i,a;2), 

where (x^jXg) is the center of the two-dimensional anti-reflection. 

More precisely, by choosing as center of the anti-reflection the flrst available data, 

along the borders, the BCs impose 

/l — ii ,12 2/1, i2 /ii + l,i2 7 /rii+ii A2 — 2/ni /rii — ii ,12 7 H 1 7 ■ ■ ■ 7 91 7 ^2 — I7 • • • 7 ^2 7 

/ii ,1 — 12 2/ii ,1 /ii ,12 + 1 7 fii ,112+12 — 2/ii ,ri2 fii ,112 — 42 7 H 1 7 ■ ■ ■ 7 '^l 7 *2 — 1 7 ■ ■ ■ 7 92- 

At the corners, the BCs impose for any ii — I, . . . ,qi and i2 ^ I, ■ ■ ■ ,q2, 



/l — ii ,1 — 12 


- 4/1,1- 


2/i,j2+i — 


2/u+i,i + /ii+i,i2+i7 


/l— 11,712+12 


= 4/1, - 


~ 2/l,„2 — 12 


2/ji_(-i,„2 + /ii + l,n2 — 12 7 


fni+ii,l — i2 


= 4/„i,i- 


~ 2/„i ,i2 + l 


2/,ii_ji,i -|- /m— 71,12+17 


/ni +11,712+12 


— 4/711,712 


2/^1, ^i_ 


-72 2/„i_ii,„2 + /7I1— 71,712— 12 



As a consequence the rectangular matrix A is reduced to a square Toeplitz-plus-Hankel 
block matrix with Toeplitz-plus-Hankel blocks, plus an additional structured low rank 
matrix. 

Moreover, under the assumption of strong symmetry of the PSF and of a mild flnite 
support condition (more precisely hi = if \ij\ > n — 2, for some j € {1,2}), the 
resulting linear system A^f = g is such that A„ belongs to the ATZ^^^ commutative 
matrix algebra [3]. This new algebra shares some properties with the r (or DST-I) 
algebra [7]- 
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Going inside the definition, a matrix An G ,47?.^^ has the following block structure 



A-n — 



Do + ^[1] 


0^ 





Di + 







Den 



t{Do,...,D,,) 















0^ 


Da + Z^^i _ 



where t(Do, ■ ■ ■ , -D^J is a block r matrix with respect to the ATZ blocks Dj^, 
ii = 1, . . . , gi and = 2 J^Hk Dt for k ^ 1, . . . ,qi. In particular, the ATl^^ block 



is associated to i\ row of the PSF, i.e., /i-^ 

as 



ID] 



2J22= — 92, 



:92 



and it is defined 





0^ 












h J- rl'^^l 
""11,92-1 T~ ^ii 

^il ,92 








^il ,92 

7j- -J- ^['^1 
"'H,<22-1 "1" '^ii 












0^ 





where z|^^ — for ^ = l,---,'?2 and T(/iij, Qj •■• i ^11,92) the unilcvcl r 

matrix associated to the one-dimensional PSF /i'^^' previously defined. 

Notice that the rank-1 correction given by the elements pertains to the contribu- 
tion of the anti-reflection centers with respect to the vertical borders, while the low 
rank correction given by the matrices pertains to the contribution of the anti- 
reflection centers with respect to the horizontal borders. 

It is evident from the above matrix structure that favorable computational properties 
are guaranteed also by virtue of the r structure. Therefore, firstly we recall the rele- 
vant properties of the two-level r algebra |3 • 

Let r„ = {An e M^(")><^("),n = (711,712), A^(n) = nin2 | A„ = Q„A„Q„} be the 
two- level T matrix algebra, i.e., the algebra of matrices that are simultaneously diag- 
onalized by the symmetric orthogonal transform 



Qn — Qni *^ Qn2i Qni — 



1 



■ sm 



stiT 



1 



(2.3) 



J s,t=l 



With the same notation as the DCT-III algebra case, the explicit structure of the 
matrix is two level Toeplitz-plus-Hankel. More precisely. 



An = Toeplitz(y) - Hankel(CT2(v^), Jcr^{V)) 
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with y = [Vb Vi ... 0...0], where each Vij, ii = l,...,qi is a the unilevel 
r matrix associated to the row of the PSF mask, i.e., Vi-^ = Toeplitz(wij) — 



Iiankel{a^{viJ,Ja^{viJ) with Vi^ = [/li^.o,- 



J , 0, . . . , 0] . Here, we denote by 
, 0] ; at the 



a'^ the double shift operator such that (T^(wiJ = [/iii,2,- 
block level the same operations are intended in block- wise sense. 
Once more, the spectral characterization is completely known since for any An € 7^ 
the related eigenvalues are given by 



(S1,S2), 



r 



where Si — l,...,rti, S2 — l,...,n2, and / is the bivariate generating function 
associated to the PSF defined in (2.2). 

As in the DCT-III case, standard operations like matrix-vector products, resolution 
of linear systems and eigenvalues evaluations can be performed by means of FST-I 
within 0{nin2 log(nin2)) (ops). For instance, it is enough to consider a FST-I applied 
to the first column of j4„ to obtain the eigenvalues 



[A„](ii,i2) = [Qn(^nei)](ii,i2)/[(3„ei](ij_i2)' *i 



,"-1,^2 = 1, 



, ,712. 



Now, with respect to the ATZ^^ matrix algebra, a complete spectral characterization 
is given in A really useful fact is the existence of a transform T„ that simulta- 

neously diagonalizes all the matrices belonging to ATZ^^, although the orthogonality 
property is partially lost. 

Theorem 2.1. Any matrix An e ATZn^. n — (ni,n2), can be diagonalized by 
Tn, i.e., 



A — T \ T 



T —T~ 



where Tn ~ Tn 



' Tn2 7 -^n 



Tn 



> T„2 , with 



T — 

rn. — 



0^ 







and T„, = 



-Qm-2P 





0^ 
0^ 



'Qm-2JP 



l-j7(m- 1), j 



. , m- 



The entries of the vector p G M™^ are defined as Pj 

J G ]gni-2xm-2 yjij^ matrix, and am is a normalizing factor chosen such that 

the Euclidean norm of the first and last column of will be equal to 1. 

Theorem 2.2. [3J Let An e ATZ^^^ , n — (711,712), the matrix related to the PSF 
hpsF — [hii,i2\ii=-qi,...,qi,i2=-q2,----q2- Then, the eigenvalues of An are given by 

• 1 with algebraic multiplicity A, 

• the 7i2 ~ 2 eigenvalues of the unilevel r matrix related to the one- dimensional 
PSF h^'^^ = Eii=-gi 'iji -<Z2: ■ ■ ■ ' Z]ji=-gi ^ji,92]' ^ac/i one with algebraic 
multiplicity 2, 

• the Til — 2 eigenvalues of the unilevel r matrix related to the one- dimensional 

jJ2T=-q ^qi^]^ each one with algebraic 



-q2 ^-qiM ; ■ 



PSF hi^} = 

multiplicity 2, 

the (ni — 2)(7i2 — 2) eigenvalues of the two-level t matrix related to the two- 
dimensional PSF hpsF- 
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Notice that the three sets of multiple eigenvalues are exactly related to the type 
of low rank correction imposed by the BCs through the centers of the anti-reflections. 
More in detail, the eigenvalues of Tn2-2{h^^^) and of Tn-^-2{h^''^) take into account 
the condensed PSF information considered along the horizontal and vertical borders 
respectively, while the eigenvalue equal to 1 takes into account the condensed infor- 
mation of the whole PSF at the four corners. 

In addition, it is worth noticing that the spectral characterization can be completely 



described in terms of the generating function associated to the PSF defined in (2.2 1 
simply by extending to the standard r evaluation grid, i.e., it holds 



m - 



1' 



where the 0— index refers to the first/last columns of the matrix 
See [21 Uj for some algorithms related to standard operations like matrix- vector prod- 
ucts, resolution of linear systems and eigenvalues evaluations with a computational 
cost of 0(riin2 log(nin2)) ops. 

It is worthwhile stressing that the computational cost of the inverse transform is com- 
parable with that of the direct transform and, at least at first sight, the very true 
penalty is the loss of orthogonality due to the first/last column of the matrix T„i. 



3. FILTERING METHODS. Owing to the ill-conditioning, the standard so- 
lution / — A^^g is not physically meaningful since it is completely corrupted by the 
noise propagation from data to solution, i.e., by the so called inverted noise. For 
this reason, restoration methods look for an approximate solution with controlled 
noise levels: widely considered regularization methods are obtained through spectral 
filtering |14l 116] . Hereafter, we consider the truncated Singular Values Decomposi- 
tions (SVDs) (or Spectral Decompositions (SDs)) and the Tikhonov (or Re- blurring) 
regularization method. 

3.1. TRUNCATED SVDs AND TRUNCATED SDs. The Singular Values 
Decomposition (SVD) highhghts a standard perspective for dealing with the inverted 
noise. More precisely, if 

is the SVD of An, i.e., C/„ and Vn are orthogonal matrices and E„ is a diagonal matrix 
with entries cri > a2 > ■ ■ ■ crj^i^n) > 0, then the solution of the linear system A^f = g 
can be written as 




where Uk and Vk denote the k column of the matrix [/„ and Vn, respectively. 
With regard to the image restoration problem, the idea is to consider a sharp filter, 
i.e., to take in the summation only the terms corresponding to singular values greater 
than a certain threshold value S, so damping the effects caused by division by the 
small singular values. Therefore, by setting the filter factors as 

, ^ f 1, if (Tfc > S, 
^'^ I 0, otherwise, 
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the filtered solution is defined as 

N{n) / T \ / T \ 

ti V ^ ti, V ^ 

Due to scale of the problem, the SVD of the matrix An is in general an expensive 
computational task (and not negligible also in the case of a separable PSF). Thus, 
an "a priori" known spectral decomposition, whenever available, can give rise to a 
valuable simplification. More precisely, let 

be a spectral decomposition of An, then the filtered solution is defined as 

where Vk and Vk denote the fc*'* column of Vn and the /c*'* row of Vn, respectively, and 
where 0^ = 1 if /c e /a, otherwise. 



3.2. TIKHONOV AND RE-BLURRING REGULARIZATIONS. In the 

classical Tikhonov regularization method, the image filtering is obtained by looking 
for the solution of the following minimization problem 

ImIi\\Anf-9\\l + ^i\\Dnf\\l (3.1) 

where /i > is the regularization parameter and D„ is a carefully chosen matrix 
(typically D„ = /„ or represents the discretization of a differential operator, properly 
adapted with respect to the chosen BCs). 

The target is to minimize the Euclidean norm of the residual ||j4„/ — 5II2 without 



explosions with respect to the quantity ||£'jia;||2. As well know, (3.11 is equivalent to 
the solution to the damped least square problem 

{AlAn + fiDlDn)f = Alg. (3.2) 

In addition, the regularization Tikhonov method can be reinterpreted in the frame- 
work of classical spectral filtering method. For instance, in the case of -D„ — /„, by 



making use of the SVD of An — ?7„S„T/^, the solution of (3.2 1 can be rewritten as 

/at = VA^-'U^g, 
where = diag((/)fe) with 0*. = CT^/(cr^ + ^), fc = 1, . . . ,iV(7i). 

A severe drawback in adopting the Tikhonov regularization approach in the case of 
An G ATZn^ is due to the fact that An ^ ATZn^, so that all the favorable compu- 
tational properties are substantially spoiled. An alternative approach, named Re- 



blurring, has been proposed in [TT1|9]: the proposal is to replace A^ by A„ in (3.2 1 



where is the blurring matrix related to the current BCs with a PSF rotated by 



180°. This approach is completely equivalent to (3.2 1 in the case of Dirichlet and 



Periodic BCs, while the novelty concerns both Reflective BCs and Anti- Reflective 
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BCs, where in general A„ ^ A^^. The authors show that the Re-blurring with anti- 
reflective BCs is computationally convenient and leads to a larger reduction of the 
ringing effects arising in classical deblurring schemes. From the modelling point of 
view, the authors motivation relies upon the fact that Re-blurring smoothes the noise 
in the right hand side of the system, in the same manner as this happens in the case 
of Dirichlet, Periodic and Reflective BCs. 

Hereafter, we consider an explanation of the observed approximation results. As pre- 
viously claimed, we focus our attention on the case of a strongly symmetric PSF, so 
that the matrix equals the matrix An. Moreover, also in this case it is evident 
that the linear system 

{Al + f,Dl)f = Ang. (3.3) 
is not equivalent to a minimization problem, again because the matrix A e ATZ^^ 



is not symmetric. Nevertheless, the symmetrization of (3.3 1 can be performed by 
diagonalization, so obtaining 

{AXn + M!),„)/ = ^A.n9, (3.4) 



where / — Tnf and g — Tng. In such a way (3.4) is again equivalent to the minimiza- 
tion problem 

mm \\AA,nTnf - TngWl + n\\AD,nTnf\\l, (3.5) 

or equivalently, again by making use of the diagonalization result, to 

mm ||f„(A„/ ~ g)\\l + ^l\\fnDnf\\l (3.6) 



Clearly, the last formulation in ( 3.6 ) is the most natural and it allows to claim that the 
Re-blurring method can be interpreted as a standard Tikhonov regularization method 
in the space transformed by means of T„. 

Recalling that r„ is not an orthogonal transformation, the goal becomes to com- 
pare ||T„/||2 and ||/||2, that is to bound ||r„||2 = ||r„J|2||r„J|2, being ||T„/||2 < 
l|r„||2||/||2. 

A quite sharp estimate of such a norm can be found by exploiting the structure of the 
unilevel matrix T™ e M"""™. Let / = [/2, . . . , fm-i], it holds that 

WTmfWl = aljf + \\Qn.-2{-fiP + f - fnJp)\\l + alf^ 

= «™(/? + /™) + II - flP +f- fnJpWl 

<ai{fi + fi) + {\\fh + m + \m\ph? 

< «™(/f + /™) + II/II2 + 3|b||^||/||^ + 4||p||2||/||^ 

<(l+2||p||2)^||/||^, 

being = 1 + ||p||2- Since, by definition, ||p||2 — m, we have 

||r™||2 < l + 2||p||2~2Vrr^. (3.7) 



Notice that the bound given in (3.7l is quite sharp, since for instance UrmeiHj equals 
l + 2Wi. 
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4. CROSS-CHANNEL BLURRING. Hereafter, we extend the analysis of 
the deblurring problem to the case of color images digitalized, for instance, according 
to the standard RGB system. Several techniques can be used for recording color im- 
ages, but the main problem concerns the fact that light from one color channel can 
end up on a pixel assigned to another color. The consequence of this phenomenon is 
called cross-channel blurring among the three channels of the image and it sums up to 
the previously analyzed blurring of each one of the three colors, named within-channel 
blurring. 

By assuming that the cross-channel blurring takes place after the within-channel blur- 
ring of the image, that it is spatially invariant and by assuming that the same within- 
channel blurring occurs in all the three color channels, the problem can be modelled 
US] as 

(^color® A„)/ =.9-?? (4.1) 

with A,, G ^N(n)-KN(n)^ (^1,713), N {v) = riina, and 



)lor 



a-rb 
agb 

dbr Clbg abb 



^gr 0,gg 



The row-entries denote the amount of within-channel blurring pertaining to each 
color channel; a normalized conservation law prescribes that Acoiore = e, e = [1 1 I]"'". 
Lastly, the vectors /, g, 7] e ]R3iV(n) ^^^.g assumed to collect the three color channels in 
the RGB order. 

Clearly, if Acoior = -^3, i-e., the blurring is only of within-channel type, the problem is 
simply decoupled into three independent gray-scale deblurring problems. 
In the general case, taking into account the tensorial structure of the whole blurring 
matrix Acoior ^ An is evident that the truncated SVDs and SDs can be formulated 
as the natural extension of those considered in the within-blurring case. Notice that 
in the case of SDs, we will consider a SVD for the matrix Acoior, since it naturally 
assures an orthogonal decomposition, no matter about the specific matrix, while its 
computational cost is negligible with respect to the scale of the problem. In addition, 
we tune the filtering strategy with respect the spectral information given only by the 
matrix An, i.e., for any fixed (or Afc) we simultaneously sum, or discard, the three 
contribution on / related to the three singular values of ^coior- 

With respect to the Tikhonov regularization method, the approach is a bit more 
involved. Under the assumption An = Aj^ = V^„A„y„, the damped least square 
problem 

[(^color ® ^n)'^(^color ® An) + nhn]f = (^color «) ^„)^5 

can be rewritten as 

[(^color^color ) ® y„A2y„ + ^(/3 ® /„)]/ = (Acolor ® VnAnVnfg. (4.2) 



Thus, by setting S^n ~ Is'^ Ki, / = S^nf, 9 — Ssng, (4.2 1 can be transformed in 

'S'3n[(ALor^color) KA,^,^ + n{h ® In)]S^n f = '53„(^color <^ VnKVnfS^^Q, 

SO obtaining the linear system 

[(^Lor^olor) ®Al + ® In)]f = (^folor ® ^n)9, 
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that can easily be decoupled into nin2 linear systems of dimension 3. 
Clearly, in the case of any matrix A„ e C„, all these manipulations can be performed 
by means of an orthogonal transformation S^n- Notice also that the computational 
cost is always 0(nin2 log nin2) ops. 

With respect to An — r„A„T„ e ATZ^^, we can consider the same strategy by 
referring to the Re-blurring regularization method. More precisely, the linear system 

[(^mlor^color) ® AI + ll{h ® /„)]/ = {A^^Xo. ® ^n)9 

can be transformed in 

[(^Lor^olor) <»Al + ® In)]f = (v^Lor ® A„).g. 

Though the transformation S^n = I?,®Tn is not orthogonal as in the Reflective case, 
the obtained restored image are fully comparable with the previous ones and the 
computational cost is still 0(nin2 log nin2)) ops. 

5. NUMERICAL TESTS. 

5.1. SOME COMPUTATIONAL ISSUES. Before analyzing the image re- 
storation results, we discuss how the methods can work without reshaping the involved 
data. In fact, the tensorial structure of the matrices, obtained by considering Reflec- 
tive and Anti-Reflective BCs, can be exploited in depth, so that the algorithms can 
deal directly, and more naturally, with the data collected in matrix form. Hereafter, 
we consider a two-index notation in the sense of the previously adopted row-wise or- 
dering. 

since Vk = vl"^' ® wL"^' is represented in 



In the SD case considered in Section 



3.1 



matrix form as the required scalar product can be computed as 



Vk9 



(~[ni]\'^ ~[n2] 



0G, 



where denotes the summation of all the involved terms after a element-wise product. 
Clearly, v^. — v^^^^ is represented in matrix form as '^["^'(^^i"^')"^- In a similar 
manner, in the case of the SVD of An with separable PSF h — hi ^ /12, we can 
represent Vk = in matrix form as and u[ — (uj."^' u^^^^)'^ 

The eigenvalues required for the SD can be stored into a matrix A* e IR"!^"^^ ^]-^g 
case of An G C„ this matrix can be evaluated as 

T 

A* 



(Vn,A*Vl ./Vn^ElViy 



where A* g u^sxni (je^otes the first column of An and El the first canonical basis 
vector, reshaped as matrices in column-wise order. In addition, the two-level direct 
and inverse transform y = VnX and y = VnX can be directly evaluated on a matrix 
data as 

Y = Vn,XVl^{Vn,{VniXff and Y =Vn,XVn,^{Vn,{Vn,Xff 

by referring to the corresponding unilevel transforms. 

In the same way, the eigenvalues required in the case of An G ATZ^ can be suitably 
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stored as 





1 


A*(r„,. 


-2m) 


1 


A* = 


A*(r„,_2(/i-)) 


A*(t„. 


-2{h)) 


A*(t„,_2(/i^)) 




1 


A*(r„,- 




1 



with reference to the notations of Theorem |2.2| where the eigenvalues of the unilevel 
and two-level r matrices are evaluated as outlined in Section 12.21 
Lastly, the linear systems obtained, for any fixed /x, in the case of Tikhonov and Re- 
blurring regularization methods can be solved with reference to the matrix <!>„ of the 
corresponding filter factors by applying the Reflective and Anti-Reflective transforms 
with a computational cost 0(nin2 log nin2) ops. 

5.2. TRUNCATED DECOMPOSITIONS. In this section we compare the 
effectiveness of truncated spectral decompositions (SDs) with respect to the standard 
truncated SVDs both in the case of Reflective and Anti-Reflective BCs. Due to scale 
of the problem, the SVD of the matrix An is in general an expensive computational 
task (and not negligible also in the case of a separable PSF). Thus, a spectral decom- 
position, whenever available as in these cases, leads to a valuable simplification. 
Firstly, we consider the case of the separable PSF caused by atmospheric turbulence 



1 



exp 



1 

' 2 I cr„ 




where and (7;, determine the width of the PSF itself. Since the Gaussian function 
decays exponentially away from its center, it is customary to truncate the values in 
the PSF mask after an assigned decay |ii|, \i2\ < I- It is evident from the quoted defi- 



nition that the Gaussian PSF satisfies the strong symmetry condition (1.7 1. Another 
example of strongly symmetric PSF is given by the PSF representing the out-of-focus 
blur 



1 

7rr^ 

0, 



if «i -f i| < 
otherwise, 



where r is the radius of the PSF. 

In the reported numerical tests, the blurred image g has been perturbed by adding a 
Gaussian noise contribution rj — rjnV with v fixed noise vector, rjn — p||g||2/||'^||2, and 
p assigned value. In such a way the Signal Noise Ratio (SNR) [,5j is given by 



SNR = 20 log 



10 



\9h 

H\2 



= 201ogiop-i {dB). 



5.2.1. GRAY-SCALE IMAGES. 



In Figure 5.1 we report the template true image (the FOV is delimited by a white 
frame), together with the blurred image with the Gaussian PSF with support 15 x 15 
and Ui-^ = (7^2 = 2 and the reference perturbation u, reshaped in matrix form. 
We consider the optimal image restoration with respect to the relative restoration 
error (RRE), i.e., ||/mt - /truc||2/||/tmo||2, where /mt is the computed approximation 
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of the true image /true by considering spectral filtering. More in detail, the RRE 
is analyzed by progressively adding a new basis element at a time, according to the 
non-decreasing order of the singular/eigen- values (the eigenvalues are ordered with 
respect to their absolute value). 

In the case of SDs (or SVDs related to a separable PSF) this can be done as described 
in Section 5.1 and, beside the preliminary cost related to the decomposition compu- 
tation, the addition of a new term has a computational cost equal to 4nin2 ops. The 
algorithm proposed in [1], that makes use of the Anti- Reflective direct and inverse 
transforms, is less expensive in the case of tests with few threshold values. 
Hereafter, the aim is to compare the truncated SVD with the truncated SD restora- 
tions both in the case of Reflective and Anti- Reflective BCs. Periodic BCs are not 
analyzed here, since Reflective and Anti-Reflective BCs give better performances with 
respect to the approximation of the image at the boundary. 

In Table |5.1| and |5.2| we report the results obtained by varying the dimension of the 
PSF support, the parameter p related to the amount of the noise perturbation and 
the variance of the considered Gaussian blur. As expected the optimal RRE worses as 
the parameter p increases and the Anti-Reflective BCs show better performances in 
the case of low noise levels. In fact, for low p values, the reduction of ringing artifacts 
is signiflcant, while the quality of the restoration for higher p values is essentially 
driven by the goal of noise flltering. Therefore, in such a case, the choice of the BCs 
becomes more an more meaningless since it is not able to influence the image restora- 
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Table 5.1 

Optimal RREs of truncated SVD and SD with reference to the true image in Figure^^^ (Gaus- 
sian blur Oi-^ = (7^2 = 2J. 









Rc 


flcctivc BCs 








PSF 




5x5 




11x11 




15x15 




21x21 


p =0 


SVD 





059164 





087402 





090742 





093856 


SD 





043754 





087400 





090746 





093867 


p =0.001 


SVD 


0.060278 





091964 


0.094468 





097034 


SD 





06027S 





091964 





094476 





097034 


p =0.01 


SVD 





091151 





11214 





11307 





11495 


SD 





091152 





11214 





11307 





11495 


p =0.05 


SVD 





11635 





13356 





13508 





13739 


SD 





11635 





13356 





13510 





13739 


p =0.1 


SVD 





13024 





14607 





14746 





15047 


SD 





13024 





14607 





14746 





15047 



Anti-Rcflcctivc BCs 



PSF 


5x5 


11x11 


15x15 


21x21 


P =0 


SVD 


0.039165 


0.064081 


0.086621 


0.087237 


SD 


0.038316 


0.063114 


0.083043 


0.083521 


p =0.001 


SVD 


0.062182 


0.094237 


0.098897 


0.10042 


SD 


0.059617 


0.089105 


0.092814 


0.094343 


p =0.01 


SVD 


0.096049 


0.12231 


0.12403 


0.12536 


SD 


0.091383 


0.11230 


0.11343 


0.11495 


p =0.05 


SVD 


0.12791 


0.15070 


0.15188 


0.15492 


SD 


0.11666 


0.13414 


0.13570 


0.13816 


p =0.1 


SVD 


0.14399 


0.16756 


0.16964 


0.17225 


SD 


0.13083 


0.14709 


0.14852 


0.15162 



Table 5.2 

Optimal RREs of truncated SVD and SD with reference to the true image in Figure \5ll\ ( Gaus- 
sian blur ai-^ = cn^ =5). 









RcflcctivG BC 










PSF 








11x11 




15x15 


21x21 


P =0 


SVD 





063387 





081274 





097351 





14634 


SD 





045365 





081274 





0963S7 





14634 


p =0.001 


SVD 





063915 





096243 





11449 





15217 


SD 





063915 





096274 





11449 





15217 


p =0.01 


SVD 





089032 





13343 





14947 





17397 


SD 





089032 





13343 





14946 





17397 


p =0.05 


SVD 





12203 





16002 





17339 





18335 


SD 





12203 





16002 





17339 





18335 


p =0.1 


SVD 





13412 





16793 





17963 





19057 


SD 





13412 





16793 





17963 





19057 



Anti-Rcflcctivc BCs 



PSF 


5x5 


11x11 


15x15 


21x21 


P =0 


SVD 
SD 


0.040214 
0.039437 


0.079543 
0.078970 


0.088224 
0.088832 


0.13686 
0.13129 


p =0.001 


SVD 
SD 


0.068197 
0.063575 


0.095808 
0.093247 


0.11522 
0.1127 


0.15767 
0.14893 


p =0.01 


SVD 
SD 


0.09412 
0.089038 


0.14482 
0.13611 


0.16825 
0.15270 


0.21148 
0.17446 


p =0.05 


SVD 
SD 


0.13553 
0.12253 


0.18563 
0.16269 


0.21006 
0.17439 


0.22962 
0.18414 


p =0.1 


SVD 
SD 


0.15010 
0.13487 


0.20164 
0.16D16 


0.22256 
0.18088 


0.23960 
0.19218 



tion quality. Some examples of restored images are reported in Figure |5.2[ 

More impressive is the fact that SDs give better, or equal, results with respect 
to those obtained by considering SVDs. This numerical evidence is really interesting 
in the case of Anti-Reflective BCs: despite the loss of the orthogonality property in 
the spectral decomposition, the restoration results are better than those obtained by 
considering SVD. Moreover, the observed trend with respect to the Reflective BCs is 
also conserved. 

A further analysis refers to the so-called Picard plots (see Figure 5.3 1, where the co- 
efficients lu^gl, or \vkg\, (black dots) are compared with the singular values CTfc, or 
the absolute values of the eigenvalues |Afe|, (red line). As expected, initially these co- 
efficients decrease faster than (Jk, or |Afc|, while afterwards they level off at a plateau 
determined by the level of the noise in the image. 

The threshold of this change of behavior is in good agreement with the optimal k 
value obtained in the numerical test by monitoring the RRE. 

Moreover, notice that the Picard plots related to the SDs are quite in agreement with 
those corresponding to SVDs. In the case of the Anti- Reflective SD we observe an in- 
creasing data dispersion with respect to the plateau, but the correspondence between 
the threshold and the chosen optimal k is still preserved. 

The computational relevance of this result is due to the significant lower computa- 
tional cost required by the Anti-Reflective SDs with respect to the corresponding 
SVDs. 

Lastly, Table [53] reports the spectral flltering results obtained in the case of Out-of- 
Focus blur by varying the dimension of the PSF support and the parameter p related 
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Fig. 5.2. Optimal restorations of trun cated SVD and SD in the case of Reflective and Anti- 
Reflective BCs with reference to Figure\5. 1\ (Gaussian blur Oi-^ = cTi^ = 2j. 
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R-SVD. k^=4653 




AR-SVD. k^=5415 















p = 0.05 



R-SD. k^=4654 



□ 5 1 1.5 2 25 




Fig. 5.3. Picard plot of truncated SVD and SD in the case of Reflective and Anti- Reflective 
BCs with reference to Figure \57l\ (Gaussian blur Oi-^ = cFi^ = 2J. 



blur). 



Table 5.3 

Optimal RREs of truncated SDs with reference to the true image in Figure \5.1\ (Out-of-Focus 



Reflective BCs 




Anti-RcflL 


ctivc BCs 




PSF 


5x5 


11x11 


15x15 


21x21 


5x5 


11x11 


15x15 


21x21 


p ^0 


0.072593 


0.084604 


0.088323 


0.096479 


0.072821 


0.085366 


0.091252 


0.099293 


p ^0.001 


0.072671 


0.085809 


0.091035 


0.10436 


0.072904 


0.086643 


0.093929 


0.10752 


p ^0.01 


0.080016 


0.12255 


0.13569 


0.15276 


0.080427 


0.12316 


0.13803 


0.15683 


p =0.05 


0.10645 


0.15365 


0.16810 


0.18777 


0.10685 


0.15571 


0.17147 


0.19172 


p =0.1 


0.12089 


0.16314 


0.17836 


0.20471 


0.12147 


0.16482 


0.17987 


0.20829 



to the noise perturbation. The RRE follows the same trend observed in the case of 
Gaussian blur. Other image restoration tests with different gray-scale images have 
been considered in pO] . 

A more interesting remark again pertains the computational cost. Since the Out- 
of-Focus PSF is not separable, but the transforms are, the use of SDs related to 
Reflective or Anti-Reflective BCs allows to exploit the tensorial nature of the corre- 
sponding transforms, both with respect to the computation of the eigenvalues and of 



18 



C. Tablino Possio 



Tn ^ iim ■■ E Ik -I M ■: I : : : - 1 1 i^n-: I ■> f : :i -.i JHu i n .1 




Fig. 5.4. True image (FOV is delimited by a white frame) and cross-channel blurred image 
with the Gaussian PSF with support 15 X 15 and = = 2 and matrix Acoior ™ j^5. . 

the eigenvectors (or of the Reflective and Anti-Reflective transforms). 



5.2.2. COLOR IMAGES IN THE CASE OF CROSS-CHANNEL BLUR- 
RING. Here, wc analyze some restoration tests in the case the template color image 
reported in Figure |5.4| by assuming the presence of a cross-channel blurring phe- 



nomenon modelled according to (4.11. The entity of this mixing effect is chosen 
according to the matrix 



)lor ^ 



0.7 0.2 0.1 

0.25 0.5 0.25 
0.15 0.1 0.75 



(5.1) 



In Figure |5.4| is also reported the cross-channel blurred image with Gaussian PSF 
with support 15 x 15 and ai^ — = 2. Notice that the entity of the cross-channel 
blurring is not negligible, since the whole image results to be darkened and the color 
intensities of the additive RGB system are substantially altered. 
Table |5.4| reports the optimal RREs of truncated SVDs and SDs obtained by vary- 
ing the dimension of the Gaussian PSF support and the parameter p related to the 
amount of the noise perturbation. It is worth stressing that we tune the filtering 
strategy with respect the spectral information given just by the matrix An, i.e., for 
any fixed Uk (or Afc) we simultaneously sum, or discard, the three contribution on 
/ related to the three singular values of Aco\or- In fact, the magnitude of singular 
values of the considered matrix Acoior does not differ enough to dramatically change 
the filtering information given just by An- Nevertheless, also the comparison with the 
restoration results obtained by considering a global ordering justifies this approach. 
The color case behaves as the gray-scale one: as expected the optimal RRE becomes 
worse as the parameter p increases and the Anti-Refiective SD shows better perfor- 
mances in the case of low noise levels. 

In addition, by referring to Figure |5.5[ we note that the truncated SVD in the case 
of Anti-Reflective BCs shows a little more 'freckles' than the corresponding truncated 
SVD in the case of Reflective BCs. Nevertheless, for low noise levels, is just the Anti- 
Reflective SD that exhibits less 'freckles' than the Reflective SD. 



5.3. TIKHONOV AND RE-BLURRING REGULARIZATIONS. By 

considering a Gaussian blurring of the true image reported in Figure |5.1| Table |5.5| 
compares the optimal RRE obtained in the case of the Tikhonov method for Reflective 




Fig. 5.5. Optimal restorations of trun cated SVD and SD in the case of Reflective and Anti- 
Reflective BCs with reference to Figure\5.4\ (Cross-channel and Gaussian Blur Oi-^ = = 2J. 
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Table 5.4 

Optimal RREs of truncated SVD and SD with reference to the true image in Figure \5^ ( Cross- 
channel and Gaussian Blur Oi-^ = cri^ = 2j. 









Reflective BC 














Ant 


-Reflective 


BC 








PSF 




5x5 


11x11 


15x15 


21x21 




PSF 


5x5 


11x11 


15x15 


21x21 


P =0 




P =0 


SVD 





078276 


0.12114 





11654 





1178 




SVD 


0.076646 


0.1006 





1098 





10646 


SD 





078276 


0.12114 





11654 





1178 




SD 


0.074953 


0.098474 





10508 





10216 


p =0.001 




p =0.001 


SVD 





078992 


0.1212 





11663 





11792 




SVD 


0.077394 


0.10639 





1111 





11002 


SD 





078992 


0.12119 





11663 





11792 




SD 


0.075727 


0.10233 





10612 





10443 


p =0.01 




p =0.01 


SVD 





10152 


0.12396 





12088 





12198 




SVD 


0.10431 


0.12695 





12624 





12779 


SD 





10152 


0.12396 





12088 





12198 




SD 


0.100S7 


0.11737 





11805 





118 


p =0.05 




p =0.05 


SVD 





12102 


0.13853 





13743 





13844 




SVD 


0.13017 


0.15075 





15063 





15166 


SD 





12102 


0.13853 





13743 





13844 




SD 


0.12127 


0.13699 





13756 





13795 


p =0.1 




p =0.1 


SVD 





13437 


0.14898 





14854 





14947 




SVD 


0.1456 


0.16516 





16626 





16647 


SD 





13437 


0.14898 





14854 





14947 




SD 


0.13507 


0.14796 





14955 





15018 



Table 5.5 

Optimal RREs of Tikhonov and Re-blurring methods and corresponding /lott with reference to 
the true image in Figure\5.1\ ( Gaussian Blur ai-^ = =2/ 



PSF 






llx 


11 


15x 


15 




21 


x21 


P =0 


R 


0.041015 


4.1O-005 


0.079044 


9e-006 


0.086386 


l.le-005 





089556 


1.6e-005 


AR 


0.034237 


l.le-005 


0.059465 


le-006 


0.078963 


le-006 


0.079805 


le-006 


p =0.001 


R 


0.050155 


0.000188 


0.087482 


5.7e-005 


0.090825 


4.3e-005 





093071 


4.9e-005 


AR 


0.048556 


0.000163 


0.085279 


4. 60-005 


0.089388 


3.3B-005 





090821 


3.3B-005 


p =0.01 


R 


0.083456 


0.005555 


0.10748 


0.001786 


0.10863 


0.001678 





11023 


0.001573 


AR 


0.083436 


0.005536 


0.10744 


0.001792 


0.10868 


0.001691 





11019 


0.001575 


p =0.05 


R 


0.12024 


0.038152 


0.12982 


0.01929 


0.13071 


0.018417 





13307 


0.017892 


AR 


0.12049 


0.038379 


0.13006 


0.01957 


0.13096 


0.018669 





1333 


0.018105 


p =0.1 


R 


0.14767 


0.06587 


0.14721 


0.039231 


0.14822 


0.038181 





15097 


0.037893 


AR 


0.14813 


0.066251 


0.14766 


0.039707 


0.14866 


0.038644 





15144 


0.038296 



BCs and of the Re-blurring method for Anti-Reflective BCs. In addition, in Table 



5.6 the same comparison refers to the case of the Out-of-Focus PSF. 
As expected, the RRE deteriorates as the dimension of the noise level or the dimen- 
sion of the PSF support increases. Notice also that the gap between the Reflective 
and Anti-Reflective BCs is reduced also for low noise levels. Further numerical tests 
can be found in [HI H] . 



Lastly, we focus our attention on the case of the color image in Figure 5.4 The 
image restorations have been obtained by considering the transformation procedure 
outlined at the end of Section |4] Despite the RREs in Table |5.7| are bigger than in 
the gray-scale case, the perception of the image restoration quality is very satisfying 
and a little less 'freckles' than in the corresponding SDs and SVDs are observed (see 



Figure 5.6). Notice, also that the lack of orthogonality in the 5*3,1 transform related 



to the Anti-reflective BCs does not deteriorate the performances of the restoration. 

6. CONCLUSIONS. In this paper we have analyzed and compared SD and 
SVD filtering methods in the case both of Reflective and Anti-Reflective BCs. Nu- 
merical evidence is given of the good performances achievable through SDs and with 
a substantially lower computational cost with respect to SVDs. In addition, the ten- 
sorial structure of the Reflective and Anti-Reflective SDs can be exploited in depth 
also in the case of not separable PSFs. 

A special mention has to be done to the fact that the loss of orthogonality of the 
Anti-Reflective transform does not seems to have any consequence on the trend of 
the image restoration results. The analysis in the case of cross-channel blurring in 
color images allows to conflrm the quoted considerations. Finally, the Re-blurring 
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Table 5.6 

Optimal RREs of Tikhonov and Re-blurring methods and corresponding fiott with reference to 
the true image in Figure\5.1\ ( Out-of-Focus blur). 



PSF 


5 




llx 


11 


15x 


15 




21x 


21 




p =0 


R 


0.031422 


0.000172 


0.05346 


6.9e-005 


0.060954 


3.5O-005 





074785 


2 


7e-005 


AR 


0.036213 


0.000302 


0.051236 


6.8e-005 


0.06683 


5.7e-005 





084482 


5 


8e-005 


p =0.001 


R 


0.034441 


0.000271 


0.061465 


0.000145 


0.073751 


0.000101 





09074 


7 


9O-005 


AR 


0.038313 


0.000402 


0.059957 


0.000138 


0.076695 


0.000126 





095274 





000106 


p =0.01 


R 


0.069647 


0.008493 


0.11361 


0.004117 


0.12881 


0.003037 





14914 





001873 


AR 


0.070384 


0.008923 


0.11404 


0.00422 


0.12982 


0.003139 





15061 





001969 


p =0.05 


R 


0.12204 


0.053687 


0.1532 


0.030719 


0.16614 


0.022121 





1876D 





01346 


AR. 


0.12256 


0.05423 


0.15402 


0.031574 


0.16739 


0.023213 





18933 





014472 


p =0.1 


R 


0.16366 


0.092379 


0.17357 


0.055919 


0.1829 


0.042944 





20323 





028803 


AR. 


0.16433 


0.093069 


0.17485 


0.057326 


0.18457 


0.044901 





20511 





031011 



Table 5.7 

Optimal RREs of Tikhonov and Re-blurring methods and corresponding fiott with reference to 
the true image in Figure\5.4\ ( Cross-channel and Gaussian Blur = = 



PSF 










ll3 


.11 




15 


xl5 




21 


(21 


p =0 


R. 


0.069148 





000203 





11508 


0.001204 





1123 


0.000717 





11335 


0.000726 


AR 


0.062854 





000102 





091232 


7O-006 





1014 


4.4C-005 





098266 


1.5O-005 


p =0.001 


R. 


0.071259 





000312 





11515 


0.001228 





11239 


0.000744 





11347 


0.000755 


AR 


0.066734 





000209 





098658 


5.SO-005 





10276 


7.7e-005 





10111 


4.5e-005 


p =0.01 


R. 


0.094871 





004975 





11896 


0.002919 





11712 


0.002421 





1182 


0.002459 


AR 


0.094458 





004841 





1144 


0.001884 





11481 


0.00184 





11507 


0.001755 


p =0.05 


R. 


0.13209 





029798 





13662 


0.015305 





13599 


0.014896 





13669 


0.014824 


AR 


0.13239 





029944 





13561 


0.014992 





13593 


0.014772 





13611 


0.014595 



R - Til li?n> - F- - -I^Miin.i 




Fig. 5.6. Optimal RREs of Tikhonov and Re-blurring methods with reference to the true image 
in Figure pl^j (Cross-channel and Gaussian blur Oi-^ = cr;, = 2 - p = 0.05J. 
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C. Tablino Possio 



regularizing method has been re-interpreted as a standard Tikhonov regularization 
method in the space transformed by means of T„. Some numerical tests highUght the 
image restoration performances, also in the case of cross-channel blurring. 
Future works will concern the analysis of effective strategies allowing to properly 
choose the optimal regularizing parameters in the Anti-Reflective BCs case. 
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